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Abstract. We prove a law of large numbers for the right boundary in the 
model of internal DLA generated by cookie random walks in dimension one. 
The proof is based on stochastic algorithms techniques. 



1. Introduction 

The model of Internal DLA (internal diffusion limited aggregation) on 1^ is a 
basic example of cluster growth model, which has the Euclidean ball as limiting 
shape, and for which we have now very good control (presumably sharp) on the 
size of the fluctuations around this limiting shape. 

It was introduced, on a mathematical level at least, by Diaconis and Fulton [8J in 
dimension one. In this case the cluster at any time is an interval, and the analysis 
of its extremities can be transposed in terms of Polya's urn, for which we know 
that the law of large numbers and the central limit theorem hold. In particular 
fluctuations around the mean value are diffusive. 

In dimension larger than or equal to 2 it has been first studied by Bramson, 
Griffeath and Lawler [6] who proved that the limiting shape is the Euclidean ball and 
that fluctuations are sublinear. Two years later Lawler [TS] showed that fluctuations 
are subdiffusive (with an upper bound in n^/^). Finally very recently Jerison, 
Levine and Sheffield [iTJ [121 [13] > and independently Asselah and Gaudilliere [U [2] 
have obtained a bound in In n in dimension 2 and in Vlnn in dimension larger than 
or equal to 3, which agrees with computer simulations. 

In its simplest form the model is defined as follows. First at time 0, the cluster 
A(0) is empty. Next if the cluster at time n is some subset A(ji) of Z*^, run a simple 
random walk and stop it when it first exit A(ji). Then A(n + 1) is the union of 
Aiji) and the site where this random walk stopped. A lot of relevant variations 
on this model can be imagined. For instance Levine and Peres [17] assume that 
random walks can start from different sources (not only the origin) . Another natural 
modification is to change the underlying graph. For instance Blachere and Brofferio 
[5] consider graphs with exponential growth and Shellef [18] looked at the case of a 
random graph, namely a supercritical cluster of percolation. 

In another direction it is also interesting to see how the analysis can be transposed 
when we change the simple random walk by some other random process. For 
instance Enriquez, Lucas and Simenhaus [10] study the case of random walks in 
random environment in dimension one and the purpose of this paper is to study the 
case of cookies random walks (also called multi-excited random walks) in dimension 
one as well. These last processes are simple examples of self-interacting random 
walks, which have been first introduced and studied by Benjamini and Wilson [4] 
and Zerner [20] . Actually our purpose is to introduce some new technique for the 
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study of Internal DLA in dimension one and then look at a specific example. The 
case of excited random walks seemed interesting to us, because it contains many 
different situations when we let the parameters of the model vary. But most of 
our analysis carry out for rather general processes, and in particular when the 
random walk is recurrent and satisfies a central limit theorem (see Theorem 12.11 
and Proposition 12. 21 below). For cookies random walks, this has been obtained by 
Dolgopyat [9 for some values of the parameters. 

Our tools are stochastic algorithms techniques which extend those used by Dia- 
conis and Fulton in their original paper. Our main result is a law of large numbers 
for the extremities of the cluster. We first present a general result in Subsection 12. II 
and then we consider in more details the case of cookies random walks in Subsection 
12.21 The proofs are postponed to the remaining sections. 



2. Model and results 

2.1. A general law of large numbers. The internal DLA process is a sequence of 
growing clusters (intervals in dimension one) which are constructed via a collection 
{X^^\k > 0)„>i, of independent random processes on Z. We assume here that 
these processes start from and only jump to the nearest neighbor, which means 
that for all n, X^"^ = 0, and for all k, xj^^ - X*."^ = ±1. We then define the 
right boundaries {dn,n > 0) of the intervals as follows. First do = 0, and for all 
n > 0, dn+i — dn + I, if hits the site (i„ + 1 before the site dn — n — 1. 

Otherwise dn+i = dn- (We also implicitly assume that our processes visit a.s. an 
infinite number of sites, which assures that dn is well defined for all n > 0.) Now 
for n > 0, set 

dn + 1 



and for a; G Z, set 



ri") = inf{fc > : Xl"^ = x}. 



Define also the function : [0, 1] — > [0, 1], by 
hn{x) = 



rp(n) rp{n) 
^[(n+2)x] ^ [(n+2)(x-l)] 



where for any u G M, [u] denotes the integer with the smallest absolute value, at 
distance strictly less than 1 from u. Note that for all n > 0, /i„ is nonincreasing, 
^n(O) = 1 and hn{l) = 0. The next result is a general law of large numbers for c?„. 

Theorem 2.1. Assume that there exists a Junction h, such that for any x e [0, 1], 

hn{x) converges to h{x), as n ^ oo. Then there exists p € [0, 1], such that Xn — > p, 
as n ^ oo, almost surely. Moreover p is characterized by the fact that h{p—) > p > 
h{p+). 

A typical situation when this theorem applies is when the X^"^^ are i.i.d. and 
transient. Indeed in this case, it is immediate that hn{x) p :— P[Ar''"' — > +oo], 
for all X £ (0, 1). Another typical situation when the theorem applies well, in the 
recurrent regime, is when there is an invariance principle: for instance if for some 
c > 0, the process defined for f > by 

y(ri) 
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converges in law in the Skorokhod space to some process {Yt,t > 0), as n — >■ oo. In 
this case, for a; S M, set 

_ r inf {t >0 : Yt>x} ii x > 
I inf{i > : < a;} if a; < 0, 

A basic observation is the following: 

Proposition 2.2. Assume that there exists c> 0, such that {Y^"'\t > 0) converges 
in law to some process {Yt,t > 0). Assume that Ti and r_i are a.s. finite. Assume 
also that h defined for all x € [0, 1], by h{x) := ¥[Tx < T^-i], is continuous. Then, 
for any x G [0, 1], hn{x) — ?> h{x), as n ^ oo. 

Note that if K is a continuous semimartingale, with a nonzero martingale part, 
then the function h defined in the proposition is automatically continuous. 

2.2. Excited random walks. We look now more precisely at the case of excited 
(also called multi-excited or cookie) random walks, which is a typical nontrivial 
example to which our general result apply, at least in most situations. However we 
stress that for some values of a parameter, Theorem 12.11 does not apply and this 
leads to interesting unsolved questions. 

We start by recalling the definition of excited random walks and we will then 
review some known results which are relevant here. 

We define an environment as the set of two sequences , « > 1) and (pi._, i > 
1) of positive reals in (0, 1). Then the law of an excited random walk (X„, n > 0) 
in this environment is defined as follows. First Xq — 0. Then, if {Tn, n > 0) is the 
natural filtration of X, we have 

r 1/2 ifx„-o 

P(X„+i =Xn + l\Tn)^l P,^+ if Xn > and L„(X„) - i 

[ Pi- if Xn < and £„(X„) = i, 

where for all x and n > 1, Ln{x) — ^{k < n : Xk — x} is the local time of the 
walk in x at time n. Moreover X always jumps to a nearest neighbor, so for all n, 
P(X„+i = X„ + 1 I = 1 - P(X„+i = X„ - 1 I Fn). This definition slightly 
differs from usual ones, where = pi- and where does not play a special 
role, but technically this does not change anything. The reason why we use this 
definition is that it gives rise to slightly richer situations. For instance the process 
can be transient only to the right. 

Now we need some assumption on the environment. We say that a sequence of 
reals > 1) is nonnegative, resp. nonpositive, if Xi > 0, resp. Xi < 0, for all 

i > 1. We will use the following hypothesis: 

Hypothesis 2.3. Each of the sequences {pi,+ — 1/2, i > 1) and {pi^- — 1/2, i > 1) 
is either nonnegative or nonpositive. 

Under this hypothesis 

oo 

a := ^(2p,,+ - 1), 
1=1 

and 

oo 

/3:=-^(2k,_-1), 

i=l 

are well defined. 
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Consider also the process F*^"^ defined for eacfi n > 1 and for all i > by 

(1) y/") := 

It was proved by Zerner [20] that under Hypothesis I2.3[ X is recurrent if, and only 
if, a < 1 and /? < 1. Note that a recurrence criterion has also been proved without 
assuming Hvpothesis 12.31 but with the other hypothesis that pi,± = 1/2 for i large 
enough (see [Hj). Later on Dolgopyat proved the following central limit thcorcm|^ 

Theorem 2.4 (Dolgopyat j9]). Assume Hypothesis \2.3\ holds. Assume also that 
a < 1 and /? < 1 . Then the process (I'/"'* ,t > 0) converges in the Skorokhod space 
toward the (a, (3) -perturbed Brownian motion (Yt,t > 0) defined by the equation: 

(2) Yt^Bt + a sup Ys+P inf K,, 

s<t «<* 

where B is a Brownian motion. 

Perturbed Brownian motions are continuous processes which were introduced by 
Le Gall and Yor [TB]. Uniqueness in law for the equation ^ was later proved by 
Perman and Werner |19| for any a < 1 and /? < 1, and pathwise uniqueness has 
been proved by Chaumont and Doney [7]. 

We note also that both Zerner's recurrence criterion and Theorem 12.41 were 
proved in the more general case when the pi^± may be random. Here we will 
only consider deterministic environments. 

If a < 1 and /3 < 1, denote by Q^"''') the law of the (a, /?)-perturbed Brownian 
motion starting from and by ^-(q,;?) the function defined for x € [0, 1] by 

where for any x, denotes the hitting time of x. Observe that hi^^ ^-^ is continuous 
and decreasing, thus has a unique fixed point p{a,j3). Since /i(q,^)(1) = = 1 — 
^(a,^)(0), we havep(a,/3) € (0, 1). Actually ft-(Q,/3)(a:) = P(Z(1 — a, 1 — 13) < l~x), 
where Z{1 — a, 1 — /3) denotes a Beta random variable with parameters 1 — a and 
1 — /3 (see Proposition 4 (iii) in [19 ). So another formula for hi^a,i3) is given by 

Jo B(l-a, 1-/3) 
with i3(-, •) the Beta function. Note that ,3)(2;) = l — /i(^ Q,)(l — a;), which implies 
that p{a, /3) = 1 — p(/3, a). Formula (|3]) shows that p{a, (3) is a continuous function 
of (a, (3) £ (— oo, 1)^, and that p{a, /3) — ?► 1, if a — ^ 1 and (3 < 1 remains fixed (resp. 
p{a, /?) — >■ 0, if /3 — >■ 1 and a < 1 remains fixed). Moreover p is monotone in (a, (3), 
in the sense that if a < a' and (3 > 13', then p{a, (3) < p{a', (3'). 
Now we consider the model of internal DLA when {X'^"\ n > 0) is sequence of i.i.d. 
excited random walks with the same law as X above. The next result gives then a 
law of large numbers for dn in different situations: 

Theorem 2.5. Assume that Hvvothesis \2.3\ holds. If [a, (3) ^ (1,1), then there 
exists p G [0, 1], such that x„ — > p almost surely. 

• // moreover a < 1 or (3 < 1, then p only depends on a and j3 and is a 
continuous function of (a, /3) in this domain. 



Actually Zerner and Dolgopyat assumed that Pi,+ = Pi,— for all i, but their proofs extend 
immediately to the present setting. 
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• If a < 1 and /? < 1, then p ~ p{a, /3) is the unique fixed point of h(^a,p), in 
particular p G (0, 1). 

• If a> 1 or (3 > I, then p = P(X„ +oo) . 

• If a = 1 and (3 < 1, then p = \ (resp. if a < 1 and (3 — 1, then p — 0). 

On the other hand, if Pi^+ — —Pi,-, for all i >\, and a = (3 = 1, then 2;„ — > 1/2 
almost surely. 

Note that the picture is almost complete. But we do not know whether x„ converges 
in general when a = /3 — 1 and whether p = P(X„ — >■ +oo) only depends on a and 
/? when a > 1 and (3 > 1. 



3. Proof of Theorem 12.11 

The proof is based on the observation that the sequence (a;„, n > 1) is solution 
of a stochastic algorithm. First, since /i„ is nonincreasing for all n > 0, /i is also 
nonincreasing. In particular it has left and right limits in all points, and there is a 
unique p £ [0, 1], such that h{p—) > p> h{p+). Next, by definition 



for all n > 1. Therefore 

Xn+1 = Xn H — H — H — for all n > 1, 

n + 6 n + 6 n + 6 

where 

(4) liT^^ti < T^ir-n-i) - 

is the increment of a martingale, and r„+i :~ hn{xn) — p. Set now ?/„ :— Xn — P- 
Then 



+ + foralln>l, 



,2 _ „,2 



where e^_|_]^ = 2e„_|_i?/„, s„_|_i = 2r„_|_i?/„ and t„ = ©(n ^). In particular the are 
still the increments of a martingale. We claim now that 

(5) limsups„+i < 0. 

n— >oo 

Indeed, let e > be fixed. First if \xn — p\ < e, then |sn+i| < 2e, since \h„\ < 1. 
Now if Xn > p + e, since hn is nonincreasing, ft,„(a;„) — p < hn{p + e) — p. But 
hn{p + e) — > + e) < p, by hypothesis. This proves that hn{xn) ^ P l£ for n 
large enough. The case Xn < p — e is similar. This proves ([5]). We deduce that 
s'n ■= max(s„, 0) converges to 0, as n — oo. Define next (z„, n > 1) by zi = yl and 



_ 2z„ ^ri+l I * 
— — T i — T i — T + ^n+l 



+t„+i for all n > 1. 



Since |ej^| < 2, it is a well known fact that z„ almost surely (see [3] or |Dul 
Theorem 3. III. 9 p. 105]). But on the other hand it is immediate by induction on n, 
that < j/^ < z„, for all n > 1. Thus ?/„—>■ almost surely, which concludes the 
proof of Theorem O □ 
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4. Proof of Proposition 12.21 
Note that Yb = a.s., so h{0) = 1 and h{l) = 0. Now fix some x e (0, 1). 
For y e R, set T{n,y) := 7'[("+2)y]/"''- Let e > be such that x + e < 1 and 
X — e > 0. Since by hypothesis Ti and T-i are a.s. finite, there exists T > 0, such 
that P[Ti > T] < e, and P[T_i > T] < e. Since F^") converges in law to Y, by 
using Skorokhod's representation theorem, we can assume that converges a.s. 
uniformly to Y on [0, T]. In other words, if 

nn,e := |sup|r/"^ -r^i <e|, 

then we can assume that V{Qn.e) — > 1, as n — ;> oo. In particular V{Qn.e) < for n 
large enough. Now, since a; + e < 1 , we have 

P[r!„,e, T{n,x)>T]<e, 

and similarly with x ~ 1 instead of x. Moreover, 

17„,, n {T{n, x) < T{n, x-l)<T}c {r,_, < T,_i_ J ■ 

Therefore, 

hn{x) =F[T{n,x) <T{n,x-l)] < 3e + ¥[nN^e,T{n,x) < T{n,x - 1) < T] 

< 3e + P[r,_, <T,_i_e] =3e + /i(x-e), 

for n large enough. We can prove similarly that hn(x) > h(x + e) — 3e, for n 
large enough. Since e can be taken arbitrarily small, and since h is assumed to be 
continuous, this concludes the proof of the proposition. □ 

5. Proof of Theorem 12.51 

• First when a <1 and /? < 1, the result follows from Theorem 12. 11 Proposition 
Oand Theorem 

• If now a > 1 or /3 > 1, then Zerner's result f20' implies that the cookie random 
walk is transient. We then immediately get that for any x £ (0, 1), hn{x) p :— 
P[Xn — !■ oo]. Thus we can apply Theorem 12.11 with h the function equal to p on 
(0, 1), and such that h{0) = 1 and h{l) = 0. 

Note that if we are in the case a > 1 and (3 > 1, then by Lemma 8 in [20] we 
have P[X„ > for all n > 0] > and P[X„ < for aU n > 0] > 0. Thus p e (0, 1). 
On the other hand if a > 1 and /3 < 1, then p — I (resp. if a < 1 and /? > 1, then 
p^O). 

• We consider next the case a — 1 and (3 < I (the case a < 1 and /3 = 1 is 
entirely similar), which belongs to the recurrent regime. We have to prove that 
a;„ — ^ 1 almost surely. For this we use the following monotonicity result due to 
Zerner: 

Lemma 5.1 ( 20] Lemma 15). Let P be the law of a cookie random walk evolving 
in an environment {pi,±, i> 1) and P' the law of a cookie random walk evolving in 
another environment {qi^±,i > 1). Assume that pi_+ > qi_+ and pi_^ > qi^-, for all 
i > 1. Then F[Ty < T^] > F'[Ty < T^], for all x < < y. 

Two cases may appear now. Either there exists M such that "Y^^Li^ Pi.+ < 1 ^ 
Sfii Pi,+ ^ for all e > 0, there exists M — A/(e) such that 1 — e < X^fii Vi,+ < 1- 
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In any case, set = Pi,-, for all i > 1, and qi,+ = 1/2, for i > AI. In addition, 
in the first case, set — pi^+, for i < M — 1 and q'^j- = pjj — e, with e < p'^j. In 
the second case, set = for i < Af. The above lemma implies that in both 
cases, the process {dn,n > 0) associated to {pi,±,i > 1) stochastically dominates 
the process {d'^,n > 0) associated to {qi,±, i > 1). Thus liminf a;„ > p{a^, f3), for all 
e > 0, where — X]f=i "^c ~^ 1) when e — >■ 0. Moreover we already saw 

in the discussion above Theorem 12.51 that p{a, [3) — >■ 1, when a — >■ 1. This proves 
that x„ converges a.s. to 1. 

• It remains to consider the case when = —pi^^, for alH > 1, and a — j3 — \. 
In this case the processes X and —X have the same law. Thus /i„(l/2) = 1/2 for 
all n > 1, and as an immediate consequence {hn{x) — 1/2) (a; — 1/2) < 0, for all 
X £ [0, 1]. In particular, we can follow the proof of Theorem 12.11 and we find that 
Sn+i ^ 0, for all n > 1. Since it was all we needed in this proof, we conclude as 
well that Xn -^1/2 almost surely. This finishes the proof of Theorem 12.51 □ 
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